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Abstract 

In this talk I first discuss the experimental evidence for multiple scat- 
tering and the properties of the underlying event. The extensive analyses 
by Rick Field of data from CDF cannot be reconciled with traditional 
wisdom conserning multiple collisions and the AGK cutting rules. Data 
seem to imply some kind of color recombination or unexpectedly strong 
effects from pomeron vertices. 

I then discuss theoretical ideas concerning the relation between mul- 
tiple collisions and unitarity: the AGK rules, IP loops, dipole cascade 
models and diffraction. 

1 Experimental overview 
1.1 Minijet cross section 

In collinear factorization the cross section for a parton-parton subcollision is 
given by 

fosubcoll j- 

dx 1 dx2f(x 1 ,p 2 L )f(x 2 ,p 2 L )- r ^-(s = xix 2 s,p\). (1) 



j ~i^,v-^x/.»v~^x, dp a 

(Note that one hard subcollision corresponds to 2 jets.) The partonic cross 
section do j dp\ behaves like l/p* for small p±, which means that a lower cutoff, 
pj_ m im is needed. The total subcollision cross section is then proportional to 
VPimmi an d f° r pp-collisions the subcollision cross section becomes equal to 
the total cross section for pj_ TO ,„ -2.5 GeV at the Tevatron and ~ 5 GcV at 
LHC [1]. Fits to data give p± m i n ~ 2 GeV at the Tevatron and slowly growing 
with energy [2]. 

In k± -factorization there is a dynamic cutoff when the momentum exchange 
k± is smaller than the virtuality of the two colliding partons, given by kj_i and 
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fc_L2 [3]. This approach gives a very similar effect. Thus we conclude that at high 
energies the subcollision cross section is much larger than the total inelastic cross 
section, which means that on average there must be several hard subcollisions 
in one event. It was also early suggested that the increase in a tot is driven by 
hard parton-parton subcollisions [4] . 



1.2 Experimental evidence for multiple collisions 
1.2.1 Multijet events 

Besides from independent multiple subcollisions, multijet events can also orig- 
inate from multiple bremsstrahlung from two colliding partons. If we study 
four-jet events the difference between these two types of events is that in a dou- 
ble parton scattering the four jets balance each other pairwise in the transverse 
momentum plane, while such a pairwise balance is not present in the multiple 
bremsstrahlung events. The Axial Field Spectrometer at the ISR proton-proton 
collider [5] studied an "imbalance parameter" 

J= 7j[(P±l+P±2) 2 + (p±3 + P± 4 ) 2 ], (2) 

and found that there is a significant enhancement of events with small values of 
J, which thus showed a clear evidence for multiple subcollisions. 

Similar, but less clear, results for four-jet events have been observed by the 
CDF [6] and DO [7] experiments at the Tevatron. A more clear signal for multiple 
collisions at the Tevatron has instead been seen in events with three jets + 7 [8]. 
Evidence for multiple collisions has also been observed in photoproduction by 
the ZEUS collaboration at HERA [9]. 



1.2.2 Underlying event 

An inportant question is whether the hard subcollisions are correlated, or if a 
high p± event just cooresponds to two jets on top of a minimum bias event. If the 
subcollisions are uncorrelated the probability, P(n) 1 for having n subcollisions 
should be described by a Poisson distribution. This implies that 

P(2) = l -P(lf. (3) 

Here the factor 1/2 is compensating for double counting. Expressed in the cross 
sections a n = P(n)a„d (where a nc i is the inelastic non-diffractive cross section) 
this gives the relation a 2 = \o\ja n< i- The experimental groups have used the 
notation 

- = \^ < 4 > 

1 a eff 

which means that cr e // = <r n d corresponds to uncorrelated subcollisions. The 
experimental results on four-jet events referred to above find, however, that a e ff 
is much smaller than a n d- Thus at ISR one finds (for jets with p± > 4 GeV) 
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(J e ff ~ 5 mb compared to a n d ~ 30 mb, CDF finds for four-jet events (p± > 25 
GeV) and 3 jets+7 the results cr e // ~ 12 mb and ~ 14 mb respectively, to be 
compared with a n d ~ 50 mb. This means that if there is one subcollision there 
is an enhanced probability to have also another one. A possible interpretation 
is that in central collisions there are many hard subcollisions, while there are 
fewer subcollisions in a peripheral collison. 

Another sign of correlations is the observation that in events with a high 
p± jet the underlying event is enhanced, the so called pedestal effect. The UA1 
collaboration at the SppS collider studied the ^-distribution in rj around a 
jet [10]. To avoid the recoiling jet they looked in 180° in azimuth on the same 
side as the jet. The result is that for jets with E± > 5 GeV the background 
level away from the jet is roughly a factor two above the level in minimum bias 
events. Similar results have been observed in resolved photoproduction by the 
HI collaboration [11]. 

1.3 CDF analysis and the Pythia model 

Rick Field has made very extensive studies of the underlying event at the Teva- 
tron (see e.g. rcf. [12]). He has here tuned the Pythia MC to fit CDF data, and 
found tunes (e.g. tune A and tune DW) which give very good fits to essentially 
all data. In particular he has looked at the E±-How, the charged particle density, 
and p^-spectra in angular regions perpendicular to a high-pj_ jet. One notice- 
able result is that the charged multiplicity in this "transverse" region grows 
rapidly with the pj_ of the trigger jet up to p±{ charged jet) w 6 GeV, and then 
levels off for higher jet energies at twice the density in minimum bias events. 
Also the charged particle spectrum has a much higher tail out to large p±_ in 
events with a high pj_ jet, compared to the distribution in minimum bias events. 
The multiple collisions have a very important effect in the MC simulations, and 
the data cannot be reproduced if they are not included. 

The version of the Pythia MC used by Field is an implementation of an 
early model by Sjostrand and van Zijl [2]. In this model it is assumed that 
high energy collisions are dominated by hard parton-parton subcollisions, and 
also minimum bias events are assumed to have at least one such subcollision. 
To be able to reproduce the observed pedestal effect, the parton distribution 
is assumed to have a more dense central region, and is described by a sum of 
two (three-dimensional) Gaussians. For fixed impact parameter, 6, the number 
of subcollisions is assumed to be given by a Poisson distribution, with an av- 
erage proportional to the overlap between the parton distributions in the two 
colliding protons. Integrated over the impact parameter this gives a distribution 
which actually can be well approximated by a geometric distribution, that is a 
distribution with much larger fluctuations than a Poisson. 

The Pythia model does not include diffraction, and describes only non- 
diffractive inelastic collisions. Diffraction is related to the fluctuations via the 
AGK cutting rules [13]. In QCD a single pomcron exchange can be repre- 
sented by a gluon ladder. The diagram for double pomeron exchange can be 
cut through zero, one and two of the exchanged pomerons, with relative weights 
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1, —4, and 2. If we add the contributions to k cut pomerons from diagrams 
with an arbitrary number of exchanged pomerons, then we get for k > 1 with 
the weights in ref [13] a Poisson distribution. For fixed impact parameter the 
assumptions in the Pythia model are thus in agreement with the AGK rules. 

1.4 Relation E± — n c h 

Although Field's tunes of the Pythia model give good fits to data, there are 
still problems. The relation between transverse energy and hadron multiplicity 
is not what has been expected. In the AGK paper a cut pomeron was expected 
to give a chain of hadrons between the remnats of the two colliding hadrons, and 
two cut pomerons should give two such chains and therefore doubled particle 
density. This is in contrast to the CDF data, where E± grows more than the 
multiplicity in multiple collision events. 

The original AGK paper was published before QCD, and based on a mul- 
tiperiphcral model. However, also in QCD the DGLAP or BFKL dynamics 
gives color-connected chains of gluons. In the hadronization process the gluon 
exchange ought to give two triplet strings (or cluster chains) stretching between 
the projectile remnants, and in the spirit of AGK two cut pomerons should give 
four such triplet strings. Field's tunes seem instead to indicate some kind of 
color recombination which reduces the effective string length. (Similar recom- 
binations have been studied by Ingelman and coworkors [14].) 

In the Pythia model used by Field different possibilities for the color con- 
nection between the partons involved are studied. The most common parton 
subcollisions are gg —>■ gg, and as mentioned above this is expected to give 
two strings between the projectile remnants. Initial state radiation gives extra 
gluons, for which the color ordering agrees with the ordering in rapidity. There- 
fore these emissions do not increase the total string length very much, and as a 
consequence they increase E±_ more than they increase the hadron multiplicity. 

From the experimental data it was noted already in ref. [2] that two sub- 
collisions could not give doubled multiplicity, as expected from four strings as 
discussed above. It was therefore assumed that the second subcollision could 
give a just single double string connecting the two outgoing gluons. Another 
option was replacing the gluons by a qq pair, connected by a single triplet string. 
This reduces the multiplicity even further. A third possibility was to assume 
that color rearrangment caused the scattered gluons to fit in the color chains 
of the first collision, in such a way that the total string length was increased as 
little as possile. This gives a minimal additional multiplicity, and in this case 
the multiple collisions have an effect on the total E± and multiplicity similar 
to the bremsstrahlung gluons (but the jets are balanced pairwise in transverse 
momentum) . The default assumption in ref. [2] was to give each of these pos- 
sibilities the same probability, 1/3. In Field's successful tunes these ratios are 
changed, such that the last option with color rcconnection is chosen in 90% of 
the cases. 

In a more recent study by Sjostrand and Skands [1] a number of improvents 
have been added to the old PYTHIA model. The hope was that with these 
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Figure 1: (a) A pomeron loop can be cut through 0, 1, or 2 of its two branches. 
This can give gaps and bumps in the particle disribution. (b) Two pomerons 
can be represented by four gluons in two color singlet pairs. Gluon exchange can 
switch the pairs (1,2) (3,4) into the singlet pairs (1,3) (2,4). A cut as indicated 
in the figure gives a localized bump in the rapidity distribution. 

modifications it would be possible to describe data without the extreme color 
reconnections which have no real theoretical motivation in QCD. Their result 
is, however, discouraging, as they were not able to tune the new model to give 
the relation beween p± and multiplicity observed in the data. 

2 Theoretical ideas 
2.1 Pomeron interactions 

We have to conclude that something important is missing in our understanding 
of high energy collisions. Although, in the AGK paper, pomeron interactions 
are assumed to give small contributions, pomeron vertices (see e.g. [15]) and 
pomeron loops may be very important. As indicated in fig. [TJi, a pomeron loop 
can give a bump in the particle density if both branches of the loop are cut, and 
a gap if the cut passes between the two branches. It is also conceivable that 
such gaps and bumps have to be included in a " renormalized" pomeron [16]. 

In QCD a pomeron is formed by two gluons in a color singlet. Two pomeron 
exchange thus corresponds to four gluons in two singlet pairs. If the pairs (1,2) 
and (3,4) form singlets, then gluon exchange can change the system so that 
instead the pairs (1,3) and (2,4) form color singlets. This corresponds to an 
effective 21P — > 2JP coupling (cf. ref. [17]). A cut with gluons 1 and 2 on one 
side and 3 and 4 on the other can then give an isolated bump in the particle 
density, as illustrated in fig.[T|3. This type of pomeron interactions can also give 
a bound state [18], which gives a pole in the angular momentum plane and an 
essential correction to the normal cut from the exchange of two uncorrelated 
pomerons. 
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Figure 2: (a) The evolution of the dipole cascade. At each step, a dipole 
can split into two new dipoles. (b) A symbolic picture of a 7*7* collision in 
y — r^-space. When two colliding dipoles interact via gluon exchange the color 
connection between the gluons is modified. The result is dipole chains stretched 
between the remnants of the colliding systems. 

We conclude that there are still many open questions. More experimental 
information is needed, and to gain insight into the dynamics it is important to 
go beyond inclusive observables, and study observables related to correlations 
and fluctuations. 

2.2 Dipole cascade models, saturation and pomeron loops 

Multiple scattering and rescattering is more easily treated in transverse coordi- 
nate space. In Mueller's dipole cascade model [19] a color dipole formed by a qq 
pair in a color singlet is split into two dipoles by gluon emission. Each of these 
dipoles can split repeatedly into a cascade, see fig. [2^,. The probability per unit 
rapidity for a split is proportional to a = N c a s /ir. When two dipole chains 
collide, gluon exchange between two dipoles implies exchange of color and a 
recoupling of the chains, as shown in fig. [2Jd. The probability for an interaction 
between two dipoles i and j, /y, is proportional to a 2 . = ir 2 a 2 /N 2 , and is thus 
formally color suppressed compared to the dipole splitting process. 

In the eikonal approximation the total scattering probability is determined 
by the expression 1 — n«(l — Aj')' which is always smaller than 1 and thus 
satisfies the constraints from unitarity. As seen in fig.[3Ji, multiple dipole-dipole 
interactions can imply that the color dipoles form closed loops, which correspond 
to the pomeron loops in fig. 

Mueller's model includes those pomeron loops, which correspond to cuts in 
the particular Lorentz frame used for the calculation, but not loops which are 
fully inside one of the colliding cascades. This implies that the formalism is not 
Lorentz frame independent, and different ways have been suggested to achieve 
a frame independent formulation (see e.g. refs. [20,21]). However, so far no 
explicitely frame independent formalism has been presented. 

In one approach the evolution is expressed in terms of interacting dipoles. 
This implies that the number of dipoles can be reduced, and the evolution of 
the projectile cascade depends on the target. Besides the 1 — > 2 dipole vertex 
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Figure 3: (a) If more than one pair of dipoles interact it can result in dipole 
loops, which correspond to pomeron loops, (b) Schematic picture of a dipole 
swing. If the two dipoles a and b have the same color, they can be replaced 
by the dipoles c and d. The result is a closed loop formed within an individual 
dipole cascade. 

there should here also be a 2 — > 1 vertex. In another approach the evolution of 
the projectile is independent of the target, and the non-interacting dipoles are 
eliminated afterwards. In this approach there is no need to reduce the number 
of dipoles in the evolution. 
Dipole swing 

A model based on the latter approach is presented in ref. [22] . In this model 
pomeron loops can be formed with the help of a recoupling of the dipole chains, 
a "dipole swing". Just as the dipole-dipole scattering, the pomeron loops in 
the cascades should be color suppressed. With a finite number of colors we can 
have not only dipoles but also higher color multipolcs. Two charges and two 
anticharges with the same color may be better approximated by two dipoles 
formed by nearby charge- ant icharge pairs. These pairs may be different from 
the initially generated dipoles, and the result is a recoupling of the dipole chain, 
as seen in fig. [3)3. The same effect can also be obtained from gluon exchange, 
which is proportional to a s and thus also color suppressed cf. to the dipole 
splitting proportional to a. 

The swing does not result in a reduction of the number of dipoles, but the 
saturation effect is obtained as the recoupled dipoles are smaller and therefore 
have smaller cross sections. Inserted in a MC the result is approximately frame 
independent, and the model describes well both the structure function in DIS 
and the pp scattering cross section [22,23], as shown in fig. [H (For these results 
also energy conservation and a running a s are very important [24].) We see here 
that the j*p cross section satisfies geometric scaling. The pp cross section is 
reduced by about a factor 4 cf. to the one pomeron exchange at the Tevatron, 
and we also see that the result of the model is the same when calculated in the 
cms as in the rest frame of the target proton, if pomeron loops are included also 
in the evolution via the dipole swing. 

Besides the total cross sections it is also possible to calculate the probability 
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Figure 4: Left: The -fp total cross section obtained from the model in The 
result is plotted as a function of the scaling variable r = Q 2 /Q s , where Q 2 S = 
Ql{x /x) x with Qo = IGcV, x = 3 • 10~ 4 , A = 0.29. Right: The total pp 
scattering cross section. Results are presented for evolution with and without 
the dipole swing mechanism. The one pomeron result and the result obtained 
in a frame where one of the protons is almost at rest are also shown. 



to have pomeron loops formed by multiple collisions in a given frame, or loops 
formed within the cascades. As examples we find at the Tevatron in the cms 
on average 2.2 loops from multiple collisions and 0.65 loops in each of the two 
cascades. In an asymmetric frame, where the total rapidity range is devided in 
4.5 + 10.5 units, we find instead 2 loops from multiple collisions, and 0.15 and 
1.35 in the two cascades respectively. In both cases this gives in total 3.5 loops. 
At LHC we obtain in the same way in total an average of 5 loops. 

Using the eikonal approximation it is besides total cross sections also possible 
to calculate elastic scattering and diffractive excitation [25], but so far it has 
not been possible to calculate exclusive final states. The aim for the future is to 
bridge the gap between dipole cascades, AGK, and traditional MC generators, 
and construct an event generator fully compatible with unitarity and the AGK 
cutting rules. 



3 Conclusions 

• Multiple collisions arc present in data. 

• Hard subcollisions are correlated. The underlying event is different from 
a minimum bias event. 

• Rick Field's tunes of the PYTHIA MC fit Tevatron data well, but the 
relation between transverse energy and multiplicity is not understood. 
This may indicate some kind of color rearrangement, or a " renormalized 
pomeron" . 
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• Multiple collisions and unitarity constraints are easier treated in trans- 
verse coordinate space. The dipole formalism can describe IP loops and 
diffraction. The application of AGK cutting rules then imply the presence 
of rapidity gaps. 

• For the future we hope to be able to combine the dipole formalism and 
traditional MC generators to obtain event generators which include diffrac- 
tion and are compatible with unitarity and AGK. 
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